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Рассмотрена задача моделирования самосогласованной динамики электронного пучка в малогабаритных бетатронах. Дается 
описание численной модели, разработанной на основе метода макрочастиц. Приведены результаты моделирования процессов 
инжекции и захвата электронов в режим ускорения в бетатронах с аксиально-симметричным и несимметричным магнитным по- 
лем. Определены оптимальные входные параметры инжекции по энергии и току пучка (20... 40 кэВ и 0,1... 1,0 А), обеспечиваю- 
щие максимальное число захваченных электронов. Численно исследованы способы повышения эффективности захвата за счет 
использования вариации внешнего магнитного поля и дополнительного отбора энергии у захваченных электронов вихревой 
тормозящей ЭДС. Это позволило увеличить коэффициент захвата с 4 до 7,4 % и захватить в ускорение до 7,4- 10 ю электронов. 


Введение 

Широкое использование малогабаритных бета- 
тронов требует существенного повышения интен- 
сивности излучения. Поэтому одной из наиболее 
важных проблем при создании бетатронов является 
повышение эффективности захвата электронов в ре- 
жим ускорения. Известно, что определяющую роль в 
механизме захвата [1,2] играет кулоновское взаимо- 
действие электронов пучка. Существующий подход 
[3, 4] исследования этого механизма захвата в клас- 
сических бетатронах, основанный на аналитических 
методах, хорошо описывает физику явления в случае 
аксиально-симметричных магнитных полей. Это да- 
ло возможность получить качественные зависимо- 
сти параметров захвата в линейном приближении от- 
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носительно малого параметра Ай/Д, в предположе- 
нии, что ширина вакуумной камеры А К значительно 
меньше радиуса равновесной орбиты Д. 

Данные допущения не применимы в малогаба- 
ритных бетатронах, в которых поперечные размеры 
пучка сравнимы с радиусом равновесной орбиты. 
Поэтому необходимо разработать методы для полу- 
чения количественных характеристик параметров 
захвата, отражающие не только влияние поля про- 
странственного заряда пучка, но и геометрию си- 
стемы, а также сложный профиль внешнего маг- 
нитного поля. Решение этой проблемы наиболее 
полно можно осуществить на основе численного 
моделирования с использованием современного 
метода макрочастиц [5]. Достоинство данного под- 
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хода заключается в получении наряду с количе- 
ственными оценками параметров захвата также де- 
тальной информации о внутренней структуре пуч- 
ка. Ниже представлена трехмерная модель захвата 
электронов в режим ускорения в цилиндрической 
системе координат {г, і, 0}. В модели используется 
полная система самосогласованной динамики 
электронного пучка, описывающая физические 
процессы в общем случае аксиально-несимметрич- 
ных магнитных полей, которые существенно влия- 
ют на эффективность захвата при многооборотной 
инжекции. Это дает наиболее адекватное описание 
нелинейной динамики пучка, что особенно важно 
для малогабаритных бетатронов. 


Описание модели 

Численная модель разработана на основе описа- 
ния электронного пучка макрочастицами [6]. Она 
построена для трехмерной области моделирования, 
совпадающей с областью ускорительной камеры 
бетатрона: 2)={-Х<?<+Г, /?,</</?,, 0<0<2л). Здесь 
К ь Я, - радиусы внутренней и внешней стенок ка- 
меры, I - аксиальная полуширина камеры. Макро- 
частицы имеют вид коаксиальных секторов (рис. 1), 
содержащих заряд 0=И е е и массу М=М е -т е , где Ы е — 
число электронов в макрочастице; е, т е - заряд и 
масса электрона. Они представляют собой укру- 
пненные частицы с удельным зарядом, равным 
удельному заряду элементарного электрона, и ха- 
рактеризуются координатами {г, і, 0} и скоростями 
{ Ѵ„ К, Ѵ д }. Размеры частиц определяются размера- 
ми расчетной сетки. Внешнее поле в общем случае 
описывается аксиальной и радиальной компонен- 
тами индукции магнитного поля В‘(г,і,Ѳ) 

и азимутальной компонентой напряженности элек- 
трического поля Е д '(г,і,Ѳ). Таким образом, модель 
дает трехмерное описание физических процессов 
как по динамике, так и по полям. 



Рис. 1. Модель макрочастицы 


Динамика макрочастиц описывается нереляти- 
вистскими уравнениями движения электрона [7]: 


сІ 2 г 
сіі 2 


— {{гѲ)В;+Е г }+гѳ 2 , 
У т о 


с! 2 Ѳ 
Б Г 



— {{Гв)в;+Е : }, 
У >Щ 


е 

У т о 



- гв;+е л 


~-(гѲ), 


( 1 ) 


где Е Ѳ Ъ =Е;+ Е ѳ - у =\Ц\-{ѴІ с) 2 -Ѵ= V 


Начальные условия определяются параметрами 
инжектора и задаются в виде: 


г=К ш , і=2 ш , Ѳ=Ѳ Ш , Ѵ=0, К=0, Ѵ Ѳ =Ѵ Ш . 
Собственное кулоновское поле пучка описывает- 
ся уравнением Пуассона для скалярного потенциала: 
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где <р=<р(г,?„Ѳ,і), р=р(г,і,Ѳ,і). 

Граничные условия удовлетворяют требова- 
ниям идеальной проводимости стенок камеры и 
периодичности по азимуту: 

<р[= 0, <р(г, 2 , 0 ) = ер{г,г,2л). 


Через потенциал вычисляются компоненты на- 
пряженности электрического поля пучка: 
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Замыкают систему уравнение для определения 
плотности заряда: 
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где М— узел сетки; О - объем ячейки сетки; л(О) - 
число макрочастиц в О. 

Моделирование самосогласованной динамики 
пучка проводится путем численного решения систе- 
мы уравнений (1-4). На каждом временном шаге по 
распределению макрочастиц в расчетной области 
вычисляется сеточная плотность заряда пучка 
р(г,і,Ѳ,(). После этого решается уравнение Пуассона 
(2) для потенциала пучка ер(г,і, Ѳ,Г), и вычисляются 
компоненты собственных полей по формулам (3). 
Уравнения движения (1) интегрируются методом 
Рунге-Кутта второго порядка точности. В соотноше- 
нии (4) за принимаем объем, порождаемый ячейкой 
сетки (рис. 1). Для решения уравнения Пуассона (2) 
используется итерационная диссипативная разност- 
ная схема [8], что позволяет сглаживать присущие 
моделям макрочастиц вычислительные высокоча- 
стотные «шумы» без искажения основных физиче- 
ских процессов. Влиянием собственного магнитно- 
го поля на движение электронов в области рассма- 
триваемых параметров можно пренебречь. 


Исследование инжекции 
и захвата электронов в ускорение 

При численном моделировании рассматривает- 
ся инжекция пучка электронов, имеющих кинети- 
ческую энергию И^ ш =20...40 кэВ. Геометрические 
параметры установки: 1=20 мм, Д,=20 мм, 
Д 2 =80 мм, і?о =45 мм. Ввод частиц моделируется ин- 
жектором, расположенным в медианной плоско- 
сти на радиусе Д,„=60 мм и имеющим аксиальную 
ширину А2 т = 1 0 мм. Компоненты внешнего маг- 
нитного поля задаются аналитически. Моделиро- 
вание проводилось для следующих конфигураций 
внешнего магнитного поля: 
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а) Аксиально-симметричное поле. Компоненты 
магнитного поля представляются в виде: 

И 2 
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где В, 0 =В,‘(К о) - значение поля В,’ на равновесной 
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ного спадания магнитного поля. 

б) Поле с азимутальной вариацией. Для четы- 
рехсекторного магнитного поля имеем: 


#” ар =В І \\ + 8В І со&{АѲ)1 В г вар = -б ,,ар — , (6) 

г 


где 8В, - глубина вариации поля. 

Результаты моделирования для 4-х вариантов 
расчетов при токе инжектируемого пучка /= 0,5 А и 
конфигурациях магнитного поля (5-6), таблица, 
представлены на рис. 2-5. На рисунках приведены 
конфигурационные портреты пучка в плоскостях 
{г, і} и {г, 9} - (а), (б), фазовый портрет пучка в 
{г, Ѵ г } - (в) и гистограммы распределения электро- 
нов в камере по энергии - (г). 


Таблица. Основные параметры численных расчетов и числа 
захваченных электронов А/ зах 


Вариант 

расчета 

Вариация 
поля 8В г , % 

Отбор энергии 
ДИ/, эВ/об. 

А/зах 'ІО'" і 

0 нс 

30 нс 

50 нс 

150 нс 

і 

0 

0 

0 

5,5 

7,1 

4,2 

2 

50 

0 

0 

8,4 

8,1 

5,1 

3 

0 

200 

0 

7,3 

9,4 

6,3 

4 

50 

200 
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9,1 

10,8 
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Рис. 2. Момент инжекции пучка, 1=5 нс 


Моделирование инжекции проводится до мо- 
мента 1=50 нс, когда накапливается предельный за- 
ряд («ІО 11 эл.), далее до /=150 нс исследуется про- 
цесс захвата в ускорение этого накопленного пуч- 
ка. При моделировании не учитывается ускорение 
пучка вихревым полем в виду малости прироста 


энергии частиц пучка, составляющую сотые доли 
процента за оборот. Соответственно, внешнее поле 
считается постоянным во времени. Это вполне 
оправданное допущение, т. к. за время установле- 
ния рассматриваемых процессов (/<200 нс), соот- 
ветствующее 50-100 об., прирост энергии частиц 
пучка пренебрежимо мал (менее 0,5 %). В процессе 
численных расчетов ведется контроль сохранения 
полной энергии системы взаимодействующих ча- 
стиц, параметры модели обеспечивали сохранение 
энергии с погрешностью менее 1 %. 


Т = 50 нс 



Рис. 3. Момент инжекции пучка, 1=50 нс 
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Рис. 4. Захват пучка при 1=100 нс 


В ходе исследований выявлен механизм захвата в 
ускорение инжектируемого пучка в бетатроне, в ко- 
тором существенное влияние оказывают поля про- 
странственного заряда пучка. Процесс захвата пучка 
в ускорение является нестационарным с выходом в 
устойчивое квазистационарное состояние, характе- 
ризуемое отсутствием потерь частиц и сохранением 
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структуры пучка в координатном и фазовом про- 
странствах (рис. 2-5). В камеру бетатрона вводится 
моноэнергетичный пучок (рис. 2, г). В процессе ра- 
диальных колебаний из-за влияния собственных ку- 
лоновских полей происходит перераспределение 
энергии между частицами пучка (рис. 3, г). Часть ча- 
стиц пучка приобретает энергию и высаживается на 
инжектор и стенки камеры, а другая часть пучка те- 
ряет энергию и захватывается в ускорение (рис. 4, г - 
5, г). Потери энергии у захваченных электронов со- 
ставляют до 15 % от энергии инжекции, это обеспе- 
чивает им радиальный отход от инжектора (рис. 5, а). 



г, см % 


Рис. 5. Захват пучка при {=150 нс 

При энергии инжекции Щ ш =40 кэВ в объеме ка- 
меры накапливается заряд, соответствующий 
Л^7-10‘°эл., захватывается в ускорение Л1 зах «4- ІО 10 эл. 
Результаты моделирования по величине захваченного 
заряда совпадают в пределах порядка с эксперимен- 
тальными данными [2, 9]. Из результатов расчетов (та- 
блица) следует, что повысить эффективность захвата 
можно путем дополнительного отбора энергии пучка 
тормозящей вихревой ЭДС и вариацией внешнего 
магнитного поля. Более существенное повышение 
эффективности захвата достигается отбором энергии 
пучка. Отбор энергии в 200 эВ/об. позволяет захватить 
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в ускорение ТѴ^б-10 10 электронов, при этом коэффи- 
циент захвата увеличивается в 1,5 раза. Совместное 
применение вариации поля и тормозящей ЭДС обес- 
печивает максимальные параметры захвата: 
7Ѵ зах «7,4'10 10 эл., ^,^«7,4 %. Коэффициент захвата рас- 
считывается по формуле [4]: К Ш =N^N Ш1І , где 7Ѵ мах - 
максимальное число электронов, удерживаемых в ка- 
мере бетатрона внешним магнитным полем. Для вы- 
шеприведенных параметров бетатрона значение 
7Ѵ мах »1 • ІО 12 электронов. 

Выводы 

1 . Разработана численная модель самосогласован- 
ной динамики электронного пучка в малогаба- 
ритных бетатронах, описывающая физические 
процессы как в аксиально-симметричных, так и 
несимметричных магнитных полях. 

2. Моделирование процессов инжекции и захвата 
в ускорение в малогабаритном бетатроне вы- 
явили оптимальные значения тока инжектируе- 
мого пучка 0,1. ..1,0 А при энергии 20. ..40 кэВ. 
При значениях тока инжекции ниже 0,1 А мож- 
но пренебречь влиянием собственных полей на 
динамику пучка, а ток 1 А является предель- 
ным, выше которого будет происходить образо- 
вание «виртуального» катода у инжектора. 

3. В ходе численных экспериментов установлено, 
что существенную роль в захвате пучка играют 
собственные кулоновские поля, обусловленные 
радиальной неоднородностью плотности заря- 
да. Получены значения предельного заряда, на- 
капливаемого при инжекции (7- 10 10 эл.) и захва- 
ченного в ускорение (4- 10 ю эл.). 

4. Исследованы возможности повышения эффек- 
тивности захвата путем вариации внешнего 
магнитного поля и отбора энергии пучка тормо- 
зящей вихревой ЭДС. Так, при энергии элек- 
тронов 40 кэВ и величине отбора энергии 
200 эВ/об., эффективность захвата повышается 
в 1,5 раза, что позволяет захватить в ускорение 
6Т0 ш эл. Максимальные параметры захвата по- 
лучены при совместном применении вариации 
внешнего магнитного поля и тормозящей ЭДС: 
А зах «7,4Т0 10 эл.,^7,4%. 
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